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PREFACE 


The  frog  muscle  samples  studied  here  were  furnished  and 
prepared  by  Professor  Alfred  Strickholm  of  the  Department  of  • 
Physiology,  Indiana  University.  He  also  supervised  their 
further  handling.  The  yeast  samples  were  prepared  by 
Dr.  Roberto  Bastos  and  Ms.  Katherine  Assimos  using  the 
facilities  of  Professor  Henry  Mahler  of  the  Chemistry  Depart- 
ment. The  E.  coli  sample  was  prepared  by  Ms.  Margaret  Loyd 
using  facilities  of  Professor  Arthur  Koch  of  the  Microbiology 
Department. 

Thanks  are  due  to  Professor  Arthur  Clouse  and 
Mr.  Robert  Addleman,  of  the  Chemistry  Department,  for  assistance 
in  setting  up  the  NMR  apparatus. 
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INTRODUCTION 


Magnetic  field  gradient  NMR  is  well  adapted  to  measuring 
diffusion  in  colloidal  systems,  including  biological  cells, 
because  the  experimental  measurement  times  by  this  method  are 
such  that  the  distances  traveled  by  the  molecules  are  of  the 
same  order  as  the  dimensions  of  the  inhomogeneities  of  the 
system.  The  result  is  that  the  apparent  diffusion  coef- 
ficients are  dependent  on  the  diffusion  time.  By  varying  this 
latter  parameter  the  dimensions  of  the  inhomogeneities  as  well 
as  the  local  diffusion  coefficients  within  them  can  be  obtained. 

The  maximum  information  is  obtained  by  the  use  of  the 
widest  possible  range  of  diffusion  times.  In  the  work  reported 
here,  a variety  of  recently  developed  technigues  have  been 
used  to  get  a much  wider  range  of  diffusion  times  than  have 
been  employed  in  previous  studies  of  diffusion  in  biological 
materials. 

The  theory  of  the  use  of  pulsed  field  gradient  NMR  to 
measure  intracellular  self  diffusion  has  been  outlined  in 
Annual  Report  No.  I1  of  this  project.  The  eguipment  necessary 
for  performing  such  experiments  has  been  assembled,  and  a 
description  was  included  in  that  report.  A preliminary 
description  of  the  results  of  a series  of  measurements  on 
various  frog  muscles  was  also  given  there. 

J.  E.  Tanner,  Self  diffusion  in  Celle  and  Tissues , NWSC/CR/ 
RDTR-6,  Naval  Weapons  Support  Center,  Crane,  IN  (June  1975). 
Avai lable  Defense  Documentation  Center  (DDC),  Cameron  Station, 
Alexandria,  V A 22314.  AD-A014601. 


In  the  intervening  time,  the  data  on  three  systems  of 
packed  single  cells  has  been  worked  up,  and  a complete  account 
will  be  given  here. 

Portions  of  the  experimental  results  have  been  presented 
at  the  March  1976  Meeting  of  the  American  Physical  Society,  at 
the  17th  Experimental  NMR  Conference  in  April  1976,  and  at  a 
symposium  on  Magnetic  Resonance  in  Colloid  and  Interface 
Science  at  the  national  meeting  of  the  American  Chemical 
Society  in  August  1976.  A preliminary  written  draft  for 
publication  in  a scientific  journal  has  also  been  completed. 

It  includes  not  only  the  work  of  the  author,  but  all  previous 
experimental  diffusion  results  from  the  literature,  using  a 
magnetic  field  gradient  method. 

During  the  course  of  the  year  the  need  has  become  clear 
for  a rigorous  derivation  of  the  relationship  between  diffusion 
coefficients,  as  measured  by  NMR,  and  diffusion  time,  for  the 
case  of  diffusion  through  periodic  barriers  of  arbitrary 
permeability.  Preliminary  work  at  deriving  such  a relationship 
will  be  indicated  here. 


THEORY 

Diffusion  in  Cellular  Systems 

The  apparent  diffusion  coefficient,  D,  of  water  in  systems 
of  biological  cells  may  be  expected  to  be  dependent  on  the 
time  allowed  for  the  measurement,  tu,  when  this  time  is  in  the 
range  10  1 to  10  seconds.  NMR  spin-echo  diffusion  measurements 
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fall  within  this  range,  and  the  expectation  of  a variable 
diffusion  coefficient  has  been  verified  a number  of  times. 

One  of  the  structures  responsible  for  the  variation  in  D 
is  the  outer  cell  membrane,  the  plasmalemma.  When  diffusion 
experiments  are  performed  in  a time  short  enough  such  that 
most  of  the  molecules  do  not  contact  this  membrane,  D repre- 
sents the  true  cytoplasmic  diffusion  coefficient,  which  we 
refer  to  here  as  D0.  As  the  time  of  the  experiment,  to,  is 
lengthened  a larger  fraction  of  the  molecules  contact  the  cell 
membranes  and  experience  their  restrictive  effects.  D decreases, 
and  as  t becomes  very  large  0 approaches  a limit,  DOT,  which 
is  readily  shown  to  be  given  by  1/D^  1 1 / D 0 + 2/aP,  where  a is 
the  average  distance  between  outer  membranes  (actually  the 
average  component  parallel  to  the  magnetic  field  gradient), 
and  P is  the  barrier  permeability.  This  relation  is  exact  in 
the  limit  that  the  membrane  thickness  is  much  less  than  a. 

It  can  readily  be  seen  that  for  impermeable  barriers 
0^=0,  whereas  for  permeable  barriers  is  finite,  but  less 

than  D . 

0 

At  intermediate  values  of  t()  it  is  obvious  that  D lies 
between  D and  D^.  Exact  mathematical  relations  for  D vs  tD, 
for  barriers  (membranes)  of  specific  geometry,  are  quite 
complicated.  Most  such  relationships  as  have  been  derived  are 
for  the  limiting  case  of  impermeable  barriers,  P=0.  There  is 
in  existence  no  rigorously  derived  relationship  valid  for 
barriers  of  arbitrary  permeability. 


Am  oxdiup  1 11  of  f ho  relationship  between  1'  a ml  t is  i|iven 


in  Figure  1.  This  is  for  the  case  of  diffusion  between 
regularly  spaced,  impermeable  planar  barriers,  where  D is 
measured  by  tne  pulsed  gradient  method,  and  the  duration  of 
the  gradient  pulses  is  small  compared  to  the  time  required  for 
the  molecules  to  diffuse  across  the  cells.  Such  is  the  case 
here,  and  this  relation  will  be  used  to  evaluate  the  data 
presented  here. 

Although  we  do  not  have  a rigorous  expression  for  the 
diffusion  coefficient  with  permeable  barriers,  we  will  make  a 
guess  so  as  to  be  able  to  treat  the  data.  We  assume  that  the 
form  of  the  expression  is  similar  to  that  for  impermeable 
barriers,  and  proceed  by  a superposition  of  the  diffusion 
coefficient  at  long  diffusion  times  and  a time-dependent 
diffusion  coefficient  as  given  for  the  case  of  impermeable 
barriers . 

Nearly  all  of  the  derivations  for  impermeable  barriers  of 
various  geometries,  with  pulsed  or  continuous  field  gradient ,• ~4 
can  be  expressed  as 

D/D0  = f(D0t,a). 

We  might  thus  write  for  permeable  barriers: 

J.  E.  Tanner  and  E.  0.  Stejskal,  J.  Chem.  Phys.  49,  1768  (1968). 
’B.  Robertson,  Phys.  Rev.  151,  273  (1966). 

‘C.  H.  Neuman,  J.  Chem.  Phys.  60,  4508  (1974), 
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log  ( 7r2D0t  /a2  ) 


Figure  1.  Relative  apparent  diffusion  coefficient  D/D0  versus  reduced 
diffusion  time,  u^D^t/a2 , for  diffusion  perpendicular  to  impermeable, 
plane  parallel  barriers.  The  two  curves  are  at  different  values  of 
y6ga.  The  upper  curve  is  appropriate  to  the  E.  col i , red  cell,  and 
yeast  measurements.  Curves  in  the  neighborhood  of  the  lower  curve 
are  appropriate  to  the  values  of  6 and  g used  for  measurements  on 
frog  muscle  at  intermediate  diffusion  times. 


i 


D/D,,  = D.„/D , + (0o-Dj/Do  x f(D  t,a). 


(1) 


A possibly  more  consistent  equation  would  be 

D/D,  = D.../D,,  ' (D  -1)  )/l)  x f[(D  -D,  ) t,a] . (2) 

Either  of  these  can  be  written  as 

D = fD0  + (l-f)D^.  (3) 

Since  f-»l  as  t-o,  and  f*o  as  t-**>  these  equations  have  the 
proper  limits  at  t=o,i>.  At  intermediate  times,  however,  they 
predict  different  behaviors. 

According  to  Equation  (1)  the  product  D t has  a fixed 

o 

value  for  a given  fractional  decrease  of  D between  D and  D^, 
independent  of  the  barrier  permeability  as  given  by  Dm.  The 
average  distance  of  diffusion  as  if  free,  x = (2D  t)1  ' , 
would  also  be  independent  of  Dn,  as  would  the  fraction  of 
molecules  which  had  encountered  the  barriers.  Thus,  (1-f) 

could  be  taken  as  representing  the  fraction  of  molecules  which 

1 /2 

had  encountered  the  barriers,  and  (D  t)  ‘ could  be  taken  as  a 

characteristic  length  for  the  system. 

According  to  Equation  (2),  the  product  (D  -0OT)t  at  a 

particular  fractional  decrease  of  D would  be  independent  of 

Dc.  Then  [(D  -D^Jt]1  ‘ could  represent  a characteristic 

length  of  the  system.  It  can  be  seen  that  x for  free 

rms 

diffusion  as  calculated  above  would  increase  as  D increased 

oo 

at  constant  f.  This  would  be  qualitatively  consistent  with 
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i 

1 tdm i 


the  assumption  that  (1-f)  represents  the  fraction  of  molecules 
which  have  been  reflected  by  the  barriers. 

I find  the  assumption  that  the  fractional  decrease  of  the 
diffusion  coefficient  depends  on  the  fraction  of  molecules 
which  have  encountered  the  barriers  more  plausible  than  that 
it  depends  on  the  fraction  reflected.  The  degree  of  reflection 
should  determine  the  magnitude  of  the  effect  on  D,  but  not  the 
time  at  which  it  sets  in.  Therefore,  Equation  (1)  will  be 
used  to  analyze  the  data  presented  here;  and  using  f as  given 
by  Figure  1.  the  barrier  distance  will  be  calculated  as 


i Wo*, 


J1 


(4) 


where  t y/z  is  the  time  when  D = (D,  + D.J/2. 

It  may  be  noted  that  Cooper,  et.  al . have  proposed 
Equation  (3)  with  the  assumption  that  (1-f)  represents  the 
fraction  of  molecules  deflected  by  the  barriers.  With  some 
inconsistency,  however,  they  then  calculate  characteristic 

4 / 2 

lengths  proportional  to  (D^t  /o) 

A Rigorous  Formula  for  Restricted  Diffusion 
with  Permeable  Barriers 

Our  uncertainty  about  how  to  analyze  diffusion  data  for 
permeable  barriers  emphasizes  the  need  for  a rigorously  derived 
expression  for  the  apparent  diffusion  coefficient  as  a function 


R.  L.  Cooper,  D.  B.  Chang,  A.  C.  Young,  C.  J.  Martin,  and 
B.  Ancker-Johnson,  Biophys.  J.  14,  161  (1974). 
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of  the  diffusion  time.  Three  approaches  for  deriving  such  an 
expression  have  been  explored  in  detail  and  at  least  two  of 
them  are  believed  to  be  feasible,  though  cumbersome. 

The  basic  problem  is  to  derive  the  distribution  function 
P(x  x,  t ) after  an  arbitrary  t for  spins  which  started  at 

O ' D 

the  same  point,  x . The  apparent  diffusion  coefficient  is 
then  calculated  from 


D = (1/t  ) In 


x. 


t ) cos[ySg(x-x  )]dxdx 


One  of  the  approaches  under  consideration  is  based  on  a 
scheme  taken  from  the  diffusion  literature  for  calculating 
P ( x | x , t ) for  a physically  appropriate  model .'  At  this 
point  it  appears  tnat  an  algebraic  expression  will  not  be 
obtained--numerical  calculations  must  be  made  for  selected 
values  of  the  parameters.  This  method  is  now  being  programmed 
for  computer  calculations. 

Another  approach  is  a modification  of  a concept  used 
earlier  to  describe  restricted  diffusion.^  Separate  algebraic 
solutions  which  are  approximately  valid  for  different  ranges 
of  the  variables  may  be  obtainable.  However,  they  are  cumber- 
some. In  addition,  there  are  problems  with  a physical  inter- 
pretation of  the  model. 


M.  Necati  Ozisik,  boundary  Value  Problems  of  beat  Conduction 

(International  Textbook  Company,  Scranton,  PA,  1963),  Chpt.  6. 

?E.  0.  Stejskal , J.  Chem.  Phys.  43,  3597  (1965). 
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The  third  approach  is  based  on  constructing  P(x  |x,  t ) 
as  a sum  of  segments  ot  a (iaussion  function.  A method  of 
performing  the  summation  has  been  found;  however,  the  validity 
of  this  approach  must  be  established. 

EXPERIMENTAL  RESULTS 

Workup  of  the  diffusion  and  relaxation  data  taken  in  the 
previous  fiscal  year  on  samples  of  human  red  cells,  E.  coli, 
and  yeast  has  been  completed.  These  results  are  reported  and 
discussed  here  along  with  the  previously  presented  results  on 
frog  muscle. 


Sample  Preparation 

The  procedure  for  preparation  of  the  samples  will  be 
briefly  summarized  here,  with  a detailed  explanation  to  be 
given  in  the  publication  presently  in  preparation. 

The  frog  muscles  were  dissected  from  freshly  pithed  frogs 
and  kept  alive  until  the  completion  of  measurements,  at  which 
time  they  were  tested  with  electrical  stimuli  and  found  to  be 
in  nearly  as  good  a condition  as  just  after  dissection.  The 
confinement  of  the  NMR  glass  sample  tube  kept  them  stretched 
to  nearly  their  full  length.  The  diffusion  measurements  were 
made  at  about  27°C,  and  in  a direction  perpendicular  to  the 
long  axis  of  the  cel  1 . 

The  human  red  cells  were  freshly  drawn  from  a healthy 
adult  into  sterile,  heparinized  tubes,  transferred  to  5 mm  NMR 
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tubes,  centrifuged  to  51  of  the  original  volume,  and  measured 
at  27 "C  over  a three  hour  period. 

The  E.  coli  culture,  strain  ML308,  was  washed  and 
centrifuged.  The  sludge  was  blotted  and  quickly  packed  into 
the  UMR  sample  tube.  The  measurements,  at  27  C,  required 
seven  hours.  The  sample  generated  about  30  of  its  volume  of 
gas  in  this  time.  The  cells  were  observed  to  be  intact  at  the 
end  of  this  period,  though  no  viability  tests  were  performed. 

Two  different  cultures  of  haploid  yeast,  strain  X2180-1B, 
were  grown.  The  first  was  centrifuged  and  blotted;  the  second 
culture  was  suction  filtered  in  several  batches  for  varying 
lengths  of  time  (degrees  of  dryness).  One  portion  was  treated 
in  a solution  of  glutaraldehyde  to  crosslink  the  protein  on 
the  cell  membrane.  A portion  of  one  sample  was  cultured  at 
the  completion  of  measurements.  The  colony  count  indicated 
1001-  viability,  though  the  uncertainty  is  perhaps  a factor  of 
two.  A similar  test  on  the  crossl inked  sample  showed  zero 
viability. 


Diffusion  Results 

The  entire  diffusion  data  from  the  frog  muscle,  E.  coli, 
and  red  cells  are  summarized  in  Figures  2 and  3,  where  the 
apparent  diffusion  coefficients  are  plotted  against  the  experi- 
mental diffusion  times.  The  apparent  diffusion  coefficients 
were  calculated  from  the  experimental  data  according  to 
In  R = -r  Dg  (A - */ 3 ) , where  R is  the  echo  attenuation  due 
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Figure  2.  Diffusion  coefficients  versus  diffusion  time  in  muscle  of 
Rana  pipiens.  Each  point  is  an  average  of  values  at  2 to  5 field 
gradient  strengths.  In  all  cases  6=0.4  msec. 

Symbols  Samples  Method 

+ 1,2  90°-180°,  alternating  sign  gradient  pulses 

X 1,2,3  90°-180°,  conventional  gradient  pulses 

Y 2,3  90°-90o-90°,  conventional  gradient  pulses 

0 4,5  90°-90°-90°,  steady  gradient 
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1 06  x D , cm2  / sec 


t 


A - 8/  3 , seconds  — •* 


Figure  3.  Diffusion  coefficient  versus  diffusion  time  in  well-packed 
samples  of  E.  coli  and  of  human  red  cells.  The  symbols  indicate 
the  type  of  pulse  sequence  used,  and  have  the  same  meaning  as  in 
Figure  2. 


to  the  gradient.  The  intensity,  duration,  and  separation  of 
the  gradient  pulses  are  given  by  g,  , and  \,  respectively.  A 
sample  of  diffusion  measurements  in  yeast  are  presented  in 
Figure  4. 

Each  of  the  points  in  the  plots  of  Figures  2 and  3 were 
obtained  from  measurements  at  two  to  five  values  of  the  field 
gradient.  For  a few  representative  measurements,  the  echo 
attenuations  are  plotted  versus  field  gradient  squared  in 
Figure  5. 

The  values  of  the  diffusion  coefficients  at  the  longest 
and  shortest  experimental  diffusion  times,  as  well  as  the  Ts 
(spin-spin)  and  T,  (spin-lattice)  relaxation  times,  are 
summarized  in  Table  1 . 

The  measurements  for  red  cells  (Figure  3)  give  the  best 
coverage  of  the  diffusion  coefficient  range.  The  diffusion 
times  extend  well  into  the  regions  of  constant  diffusion 
coefficients  Do  and  0^,  at  short  and  long  measurement  times, 
respectively.  (The  upward  trend  of  the  data  at  t >>1  sec  is 
probably  an  artifact  induced  by  the  very  low  signal  strength 
there. ) 

The  time,  t , , at  which  the  diffusion  coefficient 
reaches  its  midpoint  may  be  estimated  from  Figure  3 to  be 
6.5  msec.  We  then  calculate  a = 6.3  urn  from  Equation  (4). 

This  value  is  large,  but  still  reasonable  considering  the 
dimensions  of  the  red  cell.  We  expected  a value  larger  than 
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Figure  5.  Examples  of  the  variation  of  echo  amplitude.  A,  versus 
field  gradient  squared,  g' , at  constant  diffusion  time,  A-6/3,  for 
the  cell  samples.  The  diffusion  times  in  msec  are  indicated  on 
the  curves.  The  slopes  may  be  divided  by  3.1x1 0e  to  g.ve  the 
apparent  diffusion  coefficients.  6,  N,  g , and  x , as  well  as  A 
are  constant  for  each  curve.  ° 1 


a geometric  average  over  all  orientations  because  of  the 
deflection  effects  of  membranes  at  oblique  orientations  to  the 
field  gradient.  With  D,,  = 6.2x10”  cm  /sec  and  Du.  = 2.6x10* 
cm: /sec  we  then  calculate  P = 0.015  cm/sec  from  1/EU  = 

1/D,,  + 2/aP.  This  is  the  minimum  that  can  be  calculated  from 
these  measurements.  The  selection  of  a smaller  value  of  a 
would  yield  a larger  value  of  P. 

The  cells  of  frog  muscles  are  much  larger  than  red  blood 
cells,  and  there  was  difficulty  extending  the  diffusion  times 
to  large  enough  values  to  directly  measure  D . Assuming  that 
the  values  at  the  longest  times  are  close  to  we  have  for 
peronius  and  sartorius  muscles  0,  = 1.2xl0-5  and  l.OxlO”13 
cm' /sec,  resp.,  and  D = 1.7x10"  and  1.6x10”  cnf/sec,  resp. 
In  both  cases  we  obtain  t , * 50  msec  from  Figure  2.  Then  we 
estimate  a = 28  urn  from  Equation  (4).  This  is  reasonable, 
considering  the  known  diameter  of  about  70  um  at  the  widest 
point  in  the  cell.  We  then  calculate  P = 0.02  cm/sec.  A 
somewhat  larger  value  of  a and  smaller  value  of  D„  might  have 
been  used,  yielding  a smaller  value  of  P. 

With  semitendinosus  muscle  the  difference  between  D and 
is  smaller  than  its  precision  of  measurement,  so  that  a 
numerical  calculation  of  a and  P is  not  justified.  However, 
we  can  see  qualitatively  that  the  membrane  permeability  of 
this  muscle  is  much  larger  than  the  membrane  permeabi 1 i ties  of 
the  other  two  muscles. 
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Probably  the  major  uncertainty  in  calculating  P for  all 
of  these  muscle  samples  results  from  the  relatively  small 
change  in  D over  a very  large  change  in  diffusion  times  and 
hence  apparatus  settings.  For  the  sartorius  and  semi tendi nosus 
muscles  there  is  additional  uncertainty  due  to  a change  in 
sample  between  the  experiments  at  short  and  at  long  diffusion 
times . 

E.  coli  was  the  smallest  cell  studied,  and  here  the 
problem  was  to  obtain  measurements  at  short  enough  diffusion 
times  so  as  to  get  a good  value  for  D0.  The  cells  are  approxi- 
mately cylinders  of  diameter  1 um  and  length  3 um.  From 
Equation  (4)  we  expect  the  transition  from  D to  D^,  to  be 
centered  near  tD  = 1/2  msec.  Only  a very  slight  change  can  be 
seen  in  this  region  of  the  plot,  indicating  a very  large  water 
permeability  of  the  cell  membrane.  Taking  a = 1.0  cm, 

D,  = 5.0x10”'  cm  /sec  and  D;„  = 4.5x10”'  cm  /sec  we  calculate 
P = 1 cm/sec.  This  is  only  an  order  of  magnitude  estimate 
because  of  the  large  uncertainty  in  the  difference  D, -0m. 

Even  so  it  is  enormous.  This  measurement  should  probably  be 
repeated  with  another  sample,  dried  more  thoroughly. 

The  sharp  drop  in  diffusion  coefficient  at  longer  t, 
corresponds  to  distances  much  larger  than  the  cell  sizes.  It 
is  believed  due  to  gas  bubbles,  which  occupied  about  1/3  of 
the  sample  volume  at  the  end  of  the  experiment.  The  measurements 
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at  long  diffusion  times  were  taken  last.  An  instrumental 
artifact  is  considered  unlikely  because  the  signal  strength 
was  still  fairly  strong  here. 

Altogether,  four  different  yeast  samples  were  observed  in 
an  extensive  series  of  measurements.  The  results,  though 
unusual  in  several  ways,  were  fairly  consistent.  A typical 
example  of  the  data  is  given  in  Figure  4.  The  echo  attenuations 
are  plotted  versus  diffusion  time  times  gradient  squared.  It 
was  not  possible  to  calculate  diffusion  coefficients  from  the 
data. 

The  first  feature  which  is  striking  is  the  strong  curvature 
of  logarithm  of  echo  attenuation  versus  gradient  squared.  For 
observations  on  a single  molecular  species  in  a homogeneous 
medium  such  plots  must  be  linear.  In  inhomogeneous  media  non- 
linearities  may  arise,  but  they  are  usually  small.1'"  The 
existence  of  several  observable  species  with  markedly  different 
diffusion  coefficients  would  explain  this  particular  feature. 
However,  there  is  probably  no  molecular  species  present  in 
concentrations  comparable  to  that  of  water  in  yeast  cells. 

The  existence  of  isolated  regions  of  water  of  greatly  differing 
mobilities  is  also  not  very  plausible. 

A second  unusual  feature  of  the  results  is  that  at  any 
point  the  change  in  echo  attenuation  is  greater  for  a change 
in  diffusion  time  than  for  a change  in  gradient  squared,  i.e., 

d In  R/d  In  A >d  In  R/d  In  g‘. 
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The  manner  of  performing  the  experiments  was  to  keep  the 
relaxation  constant  by  keeping  t and  i constant  while  varying 
the  gradient,  but  to  vary  i with  t . Thus,  these  results 
could  be  mathematical ly  modeled  by  a mixture  of  two  detectable 
species  with  different  T?'s,  the  species  with  the  longer  T 
having  the  higher  diffusion  coefficient.  In  the  actual  physical 
situation,  regions  of  absorbed  and  non-absorbed  water  would 
show  just  such  a difference.  The  absorbed  water  should  have 
the  shorter  T and  also  the  lower  D.  However,  absorbed  water 
generally  constitutes  too  small  a fraction  of  the  total  ( 10») 
to  show  the  effect  observed.  Furthermore,  the  exchange  times 
between  regions  are  probably  much  shorter  than  the  experi- 
mental diffusion  times,  so  as  to  wipe  out  this  effect. 

Possible  experimental  artifacts  have  been  considered.  In 
one  run  there  was  independent  reason  to  suspect  superposition 
of  the  observed  echo  and  the  free  induction  decay  from  the 
previous  90  pulse.  After  correction  for  this,  however,  the 
unusual  features  remained. 

Thus  the  yeast  data  are  unexplained.  The  experiments 
will  be  repeated  at  the  earliest  opportunity.  A lower  temper- 
ature will  be  employed  to  insure  that  cell  degradation  is  not 
occurring  during  measurements.  Additional  apparatus  precautions 
and  cross  checks  will  also  be  employed. 

Turning  to  the  summary  of  diffusion  results  presented  in 
Table  1,  the  column  labeled  "diffusion  at  short  times"  represents 
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diffusion  coefficients  at  experimental  times  of  0.3  to  0.5  msec, 
which  is  short  enough  so  that  the  restrictive  effect  of  the 
outer  membrane  should  be  small  in  all  cases,  except  possibly 
that  of  E.  coli.  Thus,  these  values  represent  true  intra- 
cellular diffusion  coefficients. 

First  we  see  that  all  of  the  values  are  lower  than  that 
for  pure  water.  The  water  content  of  these  samples  ranges 
between  65  and  75  . Since  much  of  the  solids  consist  of  high 
molecular  weight  proteins,  cytoplasm  is  somewhat  similar  to 
ordinary  solutions  of  high  molecular  weight  polymers  in  low 
molecular  weight  solvents.  In  previous  studies  of  such 
systems'  ' J it  has  been  found  that  the  presence  of  about  30%  of 
the  soluble  polymer  decreases  the  diffusion  constant  of  a low 
molecular  weight  solvent  by  about  a factor  of  x2.  This  is 
approximately  the  effect  found  in  the  cellular  systems.  Thus 
we  feel  it  is  reasonable  to  ascribe  the  reduced  diffusion 
coefficients  within  cells  to  the  barrier  effects  of  the  high 
molecular  weight  proteins  there.  Other  investigators  (e.g. 
Abetsedarskaya , et.  a 1 . 1 ‘ ) have  argued  similarly. 

It  is  interesting  that  diffusion  is  much  more  rapid  in 
the  muscles  than  in  the  single  cells,  although  the  solids 
content  is  nearly  the  same  in  both  cases.  Apparently  the 

J.  E.  Tanner,  Macromolecules  4,  748  (1971). 

B.  D.  Boss,  E.  0.  Stejskal,  and  J.  D.  Ferry,  J.  Phys.  Chem. 

_71  , 1501  (1967). 

L.  A.  Abetsedarskaya,  F.  G.  Miftakhutdinova , and  V.  D.  Fedotov, 
Biophysics  1J3 , 750  (1968). 
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proteins  are  arranged  in  such  a manner  as  not  to  impede 
diffusion  within  muscle  cells.  This  is  plausible  considering 
the  function  of  muscle  cells. 

The  fact  that  the  values  of  D0  for  the  three  types  of 
single  cells  are  nearly  equal  is  probably  fortuitous.  Other 
single  cells  could  probably  be  found  which  would  give  signifi- 
cantly different  values. 

Relaxation  Measurements 

A number  of  measurements  of  spin-spin  (T  ) and  spin- 

2 

lattice  (T i ) relaxation  time  were  also  performed  on  the  samples. 
Some  of  these  are  presented  in  Figure  6.  The  results  are  all 
lower  than  the  value  T =T  , = 3 sec  for  pure  oxygen-free  water  at 
this  temperature.  (T  =T„  1-1/2  sec  for  air-saturated 
water.  However,  it  is  probable  that  cell  metabolism  has 
removed  most  of  the  oxygen  from  within  the  cells.)  This  is  an 
indication  of  restrictions  to  the  molecular  motion  of  the 
water.  T >>T  is  consistent  with  a small  fraction  of  molecules 

i 

with  very  restricted  motion,  perhaps  adsorbed  on  various 

inacromolecules  in  the  cells.  The  values  of  T would  be  much 

1 

less  affected  than  those  of  T by  such  a fraction  of  low 
mobility,  and  it  may  be  noted  that,  comparing  the  single  cell 
and  the  muscle  samples  with  pure  water  (Table  1 ) , T and  D are 
found  to  decrease  proportionately,  indicating  that  for  bulk 
water  molecules  the  rotational  motions  are  not  slowed  more 
than  the  diffusional  motions  by  the  other  cell  components. 
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SIGNAL,  VOLTS 


10 


t , msec 


Figure  6.  Representative  signal  decay  curves  due  to  T and  T . 

The  symbols  are:  o,  red  cells;  0,  E.  coli;  A,  yeast;  A,  frogcmuscle. 
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This  does  not  prove  the  absence  of  local  orientational  structure 
of  water,  but  shows  that  a given  water  molecule  rapidly  chanties 
its  place  in  any  such  structure  as  may  exist. 

SUMMARY 

In  summary,  we  have  measured  self  diffusion  of  water  in 
packed  single  cell  samples  of  three  different  types  and  in 
three  different  muscle  types  of  a frog.  The  measurements  were 
performed  over  a much  wider  range  of  diffusion  times  than  has 
been  done  previously.  From  the  results  at  the  shortest  times, 
the  true  intracellular  diffusion  coefficients,  uninfluenced  by 
the  outer  cell  membranes,  were  obtained.  These  diffusion 
coefficients  are  all  lower  than  that  for  pure  water,  but  by  an 
amout  which  can  probably  be  explained  by  the  obstructive 
effects  of  the  other  cellular  components. 

From  the  diffusion  coefficients  at  long  diffusion  times, 
the  permeabilities  of  the  outer  membranes  were  estimated.  The 
membrane  of  E.  coli  was  found  to  be  much  more  permeable  to 
water  than  the  membranes  of  red  cells  or  of  the  muscle  samples. 

Measurements  of  the  spin  relaxation  showed  that  the 
reorientational  motion  of  the  major  portion  of  the  cell  water 
is  not  slowed  more  than  the  translational  diffusion.  They 
also  indicate  the  presence  of  a small  fraction  of  adsorbed 
wa ter. 


Work  has  started  on  deriving  a rigorous  expression  for 
diffusion  coefficient  versus  diffusion  time  for  the  case  of 
planar  barriers  of  arbitrary  permeability. 
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Naval  Aerospace  Medical  Center 
Pensacola,  Florida  32512 


Commanding  Officer 

U.S.  Naval  Air  Development  Center 

Attention:  Aerospace  Medical  Research  Department 

Johnsville,  Warminster,  Pennsylvania  18974 

Scientific  Library 

Naval  Biomedical  Research  Laboratory 
Naval  Supply  Center 
Oakland,  California  94625 

Director,  Life  Sciences  Division 
Army  Research  Office 
3045  Columbia  Pike 
Arlington,  Virginia  22204 

Director,  Life  Sciences  Division 
Air  Force  Office  of  Scientific  Research 
1400  Wilson  Boulevard 
Arlington,  Virginia  22209 
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DISTRIBUTION  LIST  (cont.) 


ADDRESS  COPIES 

Commanding  General 

U.S.  Army  Medical  Research  & Development  Command 
Forrestal  Building 

Washington,  D.  C.  20314  1 

Professor  H.  R.  Mahler 
Department  of  Chemistry 
Indiana  University 

Bloomington,  Indiana  47401  1 

Professor  Arthur  Koch 
Department  of  Microbiology 
Indiana  University 

Bloomington,  Indiana  47401  1 

Commander 

Naval  Air  System  Command 
Department  of  the  Navy 
Washington,  D.  C.  20361 

Attention:  Code  AIR-954,  Technical  Library  2 

Code  A IR -31 0C , Dr.  H.  Rosenwasser  1 

Commander 

Naval  Sea  Systems  Command 

Naval  Sea  Systems  Command  Neadquarters 

Washington,  D.  C.  20362 

Attention:  Code  SEA-09G3,  TechnicalLibrary  2 

Code  SEA-0332,  Dr.  A.  B.  Amster  1 

Commander 

Naval  Weapons  Center 

China  Lake,  California  93555 

Attention:  Code  533,  Technical  Library  2 

Code  60,  Dr.  H.  W.  Hunter  1 

Commander 

Air  Force  Avionics  Laboratory 
Wright-Patterson  Air  Force  Base 
Ohio  45433 

Attention:  Code  AFAL/CC  1 

Commander 

Wright-Patterson  Air  Force  Base 
Ohio  45433 

Attention:  Code  AFWAL/D0,  Technical  Library  1 
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ADDRESS  COPIES 

Commander 

Naval  Surface  Weapons  Center 

White  Oak  Laboratory 

Silver  Spring,  Maryland  20910 

Attention:  Code  230,  Dr.  L.  A.  Kaplan  1 

Code  250,  Dr.  W.  McQuistion  1 

Code  WX -21 , Technical  Library  1 

Commander 

Army  Aviation  Systems  Command 
Avionics  and  Weaponization  Division 
St.  Louis,  Missouri  63166 

Attention:  Code  DRSAV-EVW  1 

Commanding  General 

U.S.  Army  Tank  Automotive  Command 

Warren,  Michigan  48090 

Attention:  Code  DRSTA-RHFL  1 

Commander 

Naval  Surface  Weapons  Center 
Dahlgren  Laboratory 
Dahlgren,  Virginia  22448 

Attention:  Code  DG-30,  Mr.  R.  Morrissette  1 

Commanding  Officer 

Frankford  Arsenal 

Philadelphia,  Pennsylvania  19173 

Attention:  Code  SARFA-MDP-Y,  Mr.  W.  Puchalski  1 

Commanding  Officer 
Edgewood  Arsenal 

Aberdeen  Proving  Ground,  Maryland  21010 

Attention:  Code  SAREA-DE-MMP,  Mr.  M.  Penn  1 

Commander 

Ballistic  Research  Laboratories 

Interrior  Ballistics  Laboratory 

Aberdeen  Proving  Ground,  Maryland  21005 

Attention:  Code  AMXBR-IB,  Mr.  J.  R.  Ward  1 

Commander 

Aeronautical  Systems  Division  (AFSC) 

Wright-Patterson  Air  Force  Base 
Ohio  45433 

Attention:  Code  ASD/ENAMC,  Mr.  M.  Edelman  1 
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ADDRESS 

Commanding  Officer 

Picatinny  Arsenal 

Dover,  New  Jersey  07801 

Attention:  Code  SARPA-FR-E-L , Mr.  T.  Boxer 

Code  SARPA-FR-E-L-C,  Dr.  F.  Taylor 
Code  SARPA-TS-S,  Technical  Library 


Commander 

Armament  Development  and  Test  Center 
Egl in  Air  Force  Base 
Florida  32542 
Attention:  Code  ADTC/SD3E 

Commander 

Air  Force  Armament  Laboratory 
Egl in  Air  Force  Base 
Florida  32542 

Attention:  Code  AFATL/DLJW,  Mr.  A.  Beach 

Commander 

Naval  Ordnance  Station 
Indian  Head,  Maryland  20640 

Commander 

Rome  Air  Development  Center 
Griffiss  Air  Force  Base 
New  York  13441 

The  Johns  Hopkins  University 

Applied  Physics  Laboratory 

■3621  Georgia  Avenue 

Silver  Spring,  Maryland  20910 

Attention:  Library  Acquisitions:  Bldg.  5,  Rm  26 

Environmental  Research  Institute  of  Michigan 
P.0.  Box  618 

Ann  Arbor,  Michigan  48107 
Attention:  IRIA  Library 

Battelle  Memorial  Institute 
TACTEC 

Columbus,  Ohio  43201 
Attention:  Ms.  Nancy  Hall 

National  Aeronautics  and  Space  Administration 
Langley  Research  Center 
Hampton,  Virginia  23665 


DISTRIBUTION  LIST  (cont.) 


ADDRESS 

National  Aeronautics  and  Space  Administration 
Lewis  Research  Center 
21000  Brookpark  Road 
Cleveland,  Ohio  44135 

Denver  Research  Institute 
Laboratories  for  Applied  Mechanics 
University  of  Denver 
Denver,  Colorado  80210 
Attention:  Mr.  Robert  M.  Blunt 

University  of  Denver 
Chemistry  Department 
Denver,  Colorado  80210 
Attention:  Dr.  John  R.  Riter 

I IT  Research  Institute 
10  W.  35th  Street 
Chicago,  Illinois  60616 
Attention:  Dr.  Elliott  Raisen 

The  Franklin  Institute 
Research  Laboratories 
Phi ladelphia , Pennsylvania  19103 
Attention:  Mr.  Gunther  Cohn 

Professor  A.  Strickholm 
Department  of  Physiology 
Indiana  University 
Bloomington,  Indiana  47401 
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